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TlCuCl3 is a quantum- spin- | system which shows a gap between the singlet ground state of the 
Cu 2+ dimers and the first excited triplet S z = +1 state for magnetic fields fioH < fioH c « 5.5 T. 
At larger magnetic fields the gap is suppressed, and a Bose-Einstein condensation (BEC) of triplets 
is supposed to occur, leading to a magnetic phase with antiferromagnetic long-range order of the 
transverse spin components. 

In this study we calculate the fraction of condensed magnetic quasiparticles of TICUCI3 from mag- 
netization M(T, if)-data. At T = K and in (i H = 6 T, this fraction is ps 98% of the total 
number of triplons. It is independent of the direction of the magnetic field and slightly decreases 
with increasing magnetic field if we assume the presence of a small intrinsic magnetic background 
with 5 = 1 magnetic moments. 

PACS numbers: 75.10.Jm, 75.30.Gw, 75.45. +j 



I. INTRODUCTION 

Low-dimensional quantum-spin systems exhibit a va- 
riety of quantum phenomena that have gained much in 
interest in the last decade 0, d, d, 0, H[ • 
TICUCI3, for example, is a material in which magnetic 
quasiparticles carrying spin 5 = 1 (spin triplet states, 
here called triplons) are believed to form a Bose-Einstein 
condensate (BEC) above a critical field hqHc ~ 5.5 T 
and at low temperatures 0, @ ■ Meanwhile several other 
materials have been found that exhibit various features 
that can be explained within the framework of a conden- 
sation of quasiparticles with integer spin [1, 0, [1, • 
The magnetic properties of TICUCI3 are determined by 
the exchange interactions between the Cu 2+ ions which 
are arranged in dimer pairs within Cu2Ci6-clusters. 
The magnetic ground state of TICUCI3 is a non-magnetic 
spin singlet that is separated from the first excited triplet 
state by an excitation gap A w 0.7 meV in zero magnetic 
field. This gap has been measured, for example, by neu- 
tron scattering and ESR measurements [10, LL1| which 
revealed that this gap is due to the strong antiferromag- 
netic interaction J= 5.68 meV in the planar dimer of 
CU2CI6. The neighboring dimers are coupled by strong 
interdimer interactions along the double chain and in the 
(1 -2) plane [13, [3. 

As soon as the external magnetic field H is larger than 
a critical field H c with g/j,BlJ.oH c (T = 0) = A (where 
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/is is the Bohr magneton and g is the Lande g- factor), 
the excitation gap closes due to the Zeeman splitting, 
and the triplet states S z = +1 are populated, eventu- 
ally forming the BEC. The 3D interdimer interactions 
drive this quantum phase transition to finite tempera- 
tures leading to a temperature dependent critical field 
H C (T). The characteristic off- diagonal long-range order 
of the BEC manifests itself in the antiferromagnetic or- 
dering of the spin system in the plane perpendicular to 
the applied magnetic field [T^. 

The idea of BEC has already been used quite successfully 
to explain the transition of " normal" to " superfiuid" 4 He 
["LH • The strong interactions that exist in liquid 4 He 
may alter the nature of the transition, however. For in- 
stance, while 90-95% of the particles of an atomic en- 
semble are in the "superfiuid" phase below the transition 
temperature of an atomic BEC, just a few percent (~ 9%) 
of the Helium-atoms are condensed in superfiuid 4 He. 
In this paper we focus on the condensed phase of triplons 
in TICUCI3 at magnetic fields n$H c < (IqH < 9 T and 
at temperatures down to T= 1.9 K. From magnetiza- 
tion M(T, H) measurements we extract the density of 
condensed triplons at T = K for H || b and H || [201]. 
Taking various possible contributions to the total magne- 
tization into account, we show that the density of triplons 
forming the condensate is in fact the same for both di- 
rections @ . We also determine the magnetic- field depen- 
dence of the fraction of triplons forming the condensate. 
The quantitative results presented here confirm the sce- 
nario of the formation of a weakly interacting Bose gas 
of triplons right above H c [5j, [l6| , and we conclude that 
the interaction increases with increasing particle density, 
i.e. with increasing magnetic field H. 
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FIG. 1: The magnetic susceptibility x(T) of TlCuCl 3 with 
HoH = 1 T applied along the crystallographic &-axis and along 
[201]. The corresponding fits to the data are discussed in the 
text 



II. THE MAGNETIC SUSCEPTIBILITY FOR 

T > 20K 

Magnetic-susceptibility measurements were performed 
in a commercial PPMS (Physical Property Measurements 
System, Quantum Design) on a TICUCI3 single crystal 
with mass m = 12.36 mg, for 1.9 K < T < 300 K at fi H 
= 1 T for H || b and H || [201]. 

The susceptibility x{T) of TICUCI3 is typical for a low- 
dimensional spin gap system, showing a well pronounced 
maximum at T Xmax ~ 36 K and an exponential de- 
crease at low temperatures indicating the existence of 
a gap A between the ground state and the first excited 
triplet state, see Fig. [T] For Heisenberg spin systems 
with identical spin subsystems that are weakly coupled 
to each other, a good fit to the data in the paramagnetic 
regime is provided by the molecular mean-field theory 
(MFT) and its extensions [l7|. We therefore used this 
approach within the model of dimers coupled by an ef- 
fective interdimer coupling J, representing the sum over 
all exchange coupling constants Jki for a given dimcr 



k interacting with neighboring dimers I [17j. An ad- 
ditional temperature-independent diamagnetic term xo 
which contains the orbital diamagnetic core contribution 
^.core (j nc i uc ii n g the background contribution of the sam- 
ple holder) and a paramagnetic Van Vleck contribution 
X VV is also considered. 

Thus we fitted the magnetic susceptibility x(T) for tem- 
peratures T > 20 K according to 



X{T) = Xo + Xmf(T), 



with 
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is the susceptibility of a non-interacting spin-dimer sys- 
tem with single spins S = h and the intradimer coupling 
J- Xmf accounts for the mean-field correction [l7| . 
For a given measured data set of x(T)-data we there- 
fore used four fitting parameters: xo> 9, J and J. We 
forced the values for xo, J and J to be identical for 
both magnetic-field directions. This restriction is physi- 
cally reasonable, since these three fitting parameters are 
independent of the magnetic-field orientation. A small 
anisotropy of the g-factor was considered, however, al- 
though it is not predicted by ESR measurements pdj . 
In fact, the obtained g-values for the two investigated 
crystallographic directions are the same within the error 
margin, see Table [U The best obtained fits are shown in 
Fig. [T] They yield a good description of the experimental 
data for T > 20 K. However, the distinct upturn in xCO 
at lower temperatures is not at all reproduced by the 
fits. We believe that this term is intrinsic for TICUCI3 
[Tl| and we shall discuss it in more detail in Chapter IIIII 
Note that the inclusion of a Curie-like term for fitting 
the data at T > 20 K does not significantly change the 
results presented in Table fl] 



TABLE I: The extracted fitting parameters from x{T) data of 
TlCuCl 3 (T > 20 K). xo, J and J were forced to be identical 
for the two field directions. 
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The value of the intradimer coupling J is close to to 
the result obtained by neutron scattering measurements 
(J/fcs - -64 K [HI). 

Unfortunately, the fit does not allow us to distinguish be- 
tween the individual interaction coupling constants con- 
tributing to J in TICUCI3 as defined in Ref.[I3], but the 
fact that J « J clearly shows the strong 3D coupling be- 
tween the dimers. From Eq. (Ila|) we find that the peak 
value Xmax = x{T ~ 36 K) increases with either increas- 
ing the g-factor or the intradimer coupling constant J, or 
by decreasing the interdimer coupling constant J. Since 
our value of J is consistent with published data from neu- 
tron scattering measurements [lOj, the slight overshoot 



3 



of the fitting curve with respect to the measured data 
around Xmax implicates an underestimate of J and/or 
an overestimate of the Lande g-values, respectively. The 
latter scenario is supported by comparing our results to 
high precision ESR- measurements [Tlj which obtain a 
value of g = 2.06 for both magnetic-field directions. 



III. THE LOW-TEMPERATURE 
MAGNETIZATION 

In the theory for a BEC of magnetic quasiparticles 
in insulating materials, the total magnetization (to be 
more precise, the total magnetic moment) M = g/J-sN is 
proportional to the total number of excited triplons N, 
which depends on both the temperature T and magnetic 
field H [5| . We therefore decided to analyze in detail the 
low-temperature region of the magnetization for both low 
magnetic fields (1 T < [IqH < [IqH c ) and high magnetic 
fields (fioH c < (j,qH < 9 T) using a consistent approach 
including adequate contributions for the respective mag- 
netic field regions. We note here that all the magnetiza- 
tion data presented in this work are expressed as mag- 
netic moment M per single Cu 2+ ion . The later used 
quantity m(T) = M(T)/N d = g^ B n(T) (where N d is 
the number of dimers and n(T) — N(T)/Nd is the total 
triplon density) differs from that by a factor 2. All values 
extracted from fits and calculations are presented in the 
latter units. 



A. The magnetization M(T) for IT < /i H < fx H c 

Fig. [2] shows the variation of M of TICUCI3 at low 
temperatures along the crystallographic b- axis and the 
[201]- direction, respectively, for magnetic fields up to 
HoH = 5 T. The magnetization decreases exponentially 
to almost zero with decreasing temperature for both crys- 
tallographic directions, but showing an upturn at low 
temperatures for low magnetic fields. With increasing 
magnetic field the anisotropic behavior of the magneti- 
zation in the two different field orientations becomes ap- 
parent. Because for both field directions the upturn in 
M at low temperatures is gradually suppressed with in- 
creasing H, the magnetization curves for H || [201] cross 
at T cross ~ 3.2 K . For H |j b a similar crossing of M(T) 
data cannot be seen in the analyzed temperature range, 
but by extrapolating the respective magnetization curves 
to lower temperatures a T cross below 2 K seems to be 
plausible. 

This crossing of M(T)-data is caused by the fact that 
the upturn in M(T) at low temperatures does not grow 
linearly with H. Moreover, this upturn is quantita- 
tively different in the two considered field directions. If 
this low-temperature contribution was due to an extrin- 
sic paramagnetic impurity phase it could be expected 
to be isotropic. We therefore consider this behavior to 
be intrinsic to TICUCI3. This upturn in M{T) can be 
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FIG. 2: The magnetization M{T) of TIC11CI3 for IT < fj, H < 
5 T applied along the crystallographic 6-axis (left) and along 
[201] (right). 



expressed as a temperature and field-dependent Curie- 
Weiss-like term that is proportional to the Brillouin func- 
tion Bs(x) with x — g/iB^oHS/ksT and a constant 
Cs- Because we assume this term to be intrinsic to the 
here studied dimer-system, it is reasonable to assign it 
to magnetic moments associated with the triplet states 
with S — 1. A similar observation confirming this fact 
was reported in [ll[. The magnetic- field dependence of 
our low-temperature data showing an almost saturated 
behavior in M(H) for H — > H c> see Fig. [3l can indeed 
be qualitatively well reproduced by a magnetization term 
that is proportional to a Brillouin function B$ (x). The 
additional ii"-dependence as observed for H < H c and 
H || [201] can be explained by taking again a diamagnetic 
term mdia — Xa ■ H and an additional paramagnetic term 
rn HL (H) (to be discussed below) into account. However, 
the quality of corresponding fits to our low-temperature 
M (Zf)-data does not allow us to clearly distinguish be- 
tween S=l and S = i. Therefore, we will consider in 
the following both scenarios for the Curie- Weiss-like term 
Cs ■ Bs(x), and we will later argue that only the 5 = 1 
case fits to our data in a physically meaningful way. 

From the expression for the free energy per unit length 
of a Heisenberg ladder [l8| 



1 + 2cosh ( MB*>HS\ 



with 



z(T) 



1 

2^ 



e k B T dk, 



z{T) (2a) 



(2b) 



we can estimate the magnetization per dimer by multi- 
plying / from Eq. (|2"a|) with a characteristic length a. 
This quantity has been suggested to correspond to an 
average lattice constant a — (a ■ b ■ c ■ sin/3) < - 1 / 3 - ) = 0.79 
nm [191 ]. where a, b and c and (3 = 96.32° are taken from 
crystallographic data of TICUCI3 [l3| . 
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FIG. 3: The magnetization M(H) of TlCuCls for T = 1.9 
K with H applied along the crystallographic 6-axis (left) and 
along [201] (right). 



Using a simple quadratic approximation for the triplon 
dispersion relation, oc A + ft 2 k 2 / '2m* (where m* corre- 
sponds to the effective mass of the triplons) , one obtains 



z(T) 



1 



e k b 1 



2^F \2m*k B T / 
For the magnetization per dimer we therefore have 
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FIG. 4: The magnetization M(T) of TlCuCls for IT < fi H < 
5 T applied along the crystallographic 6-axis (upper frames) 
and along [201] (lower frames) . The continuous and the 
dashed lines denote fits for S = 1 and S = |, respectively 
(see text). 



with 
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In order to analyze the upturn in M(T) at low T, we 
include the above mentioned magnetization term 



m , 



,(T) = 9f x B S -Cs-Bs 



V k B T 
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for fixed magnetic field H. We distinguish here between 
S = s (non-intrinsic paramagnetic impurities) and the 
scenario S — 1 (intrinsic term related to triplet states). 
We fitted the magnetization data at low enough temper- 
atures (T < 5 K) and 1 T< fi H <5 T for both field 
directions according to 
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with g = 2.06. Because we were using the T-independent 
diamagnetic contribution rridia — Xo ' H extracted from 
the high-temperature susceptibility fits presented above, 
only the gap A, the constant Cs (for S'=l or |) and 
the effective mass of a triplon m* were fitting parame- 
ters. The corresponding fits to the magnetization data 
are shown in Fig. 31 while the corresponding results for 
the fitting parameters are presented in Fig. [5J The values 
for the gap A slightly vary with magnetic field for both 
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FIG. 5: The parameters A, m* and Cs extracted from fits 
to the low-temperature magnetization data for H | 6-axis 
(upper frames) and for H || [201] (lower frames). Different 
scenarios for 

TTliip Eire indicated with open marks (.5=1) and 
filled marks (S = i), respectively. The error bars of Cs are 
of the order of the dot size. 



field directions around A/ks ~ 13 K, which is somewhat 
larger than A w 0.7 meV = 8.3 K (TTJ determined by 
neutron scattering. The triplon mass to* ps 0.2 • 10~ 29 
kg is an order of magnitude smaller than compared to 
the results from calculations and a corresponding analy- 
sis of high- field magnetization data within the Hartree- 
Fock approximation [l!|. This discrepancy might be ex- 
plained by our choice of a or by the use of the simplified 
quadratic energy-dispersion relation for this temperature 
region. The range of validity of a quadratic approxima- 
tion is indeed restricted to lower temperatures (T <1 K) 
[IS H HJ that are not accessible in our experiment. 

The Curie-like contribution Cs decreases for both cases 
5=1 and 5 = | with increasing magnetic field H. For 
both field directions Cs(H) shows a similar trend, al- 
though the variation with H is much less pronounced for 
5=1. Since Cs is expected to be a constant for a given 
magnetic-field direction, this fact already here strongly 
supports an 5 = 1 scenario for a correct description of 
the paramagnetic background. We may speculate that 



this Curie-like term with 5= 1 comes from a contribution 
of defects in the crystal or from dimers that are situated 
near the crystal boundaries. 



B. The magnetization M(T) for 5.5T < /i H < 9T 

The temperature dependence of the magnetization 
M(T) along the applied magnetic field H shows a cusp- 
like minimum at a critical temperature T C (H) for fixed 
magnetic fields H > H c , see Fig. [5] The increase of 
M for T < T c is a consequence of the condensation of 
the magnetic quasiparticles and the increasing number 
of particles N c in the ground state forming the conden- 
sate. Theoretical arguments suggest within a simplified 

model a T-dependence of M oc (l - Z-\ 2 for T < T c @ 

which is not observed in the experimental data, however. 
At high magnetic fields we have therefore fitted the low- 
temperature magnetization per dimer m(T) according 
to a more general power-law including the diamagnetic 
contribution that we extracted from high-temperature 
magnetic-susceptibility measurements, and again a net 
paramagnetic moment m up (T, H) assumed to be propor- 
tional to the Brillouin function Bs{x) for 5=1 and 
5 = i, respectively. For 6 T< /iqH <9 T we use 



M(T) 
N d 



9VB- 



N(T) 



mdt, 



no 1 - =- 



(7) 

where n cr n — N(T — T c )/Nd is the critical density at 
which condensation occurs, corresponding to the normal- 
ized magnetization m(T = T c ) = g[iBn cr it- The physical 
meaning of the exponent a, see Fig. [TJb), is not further 
discussed here, see Appendix. 

At zero temperature we have for fixed magnetic field H 



m(T = 0) = gnB(n crit + n ) + m up (T = 0) + m dl , 



(8) 



= g^i B n{Q) + m up (T = 0) + m dla 

with the total triplon density at T — 0, n(0) = n C rit +?^o 

For an ideal Bose gas n(0) corresponds to the condensate 
density n c (0). As soon as interactions between the par- 
ticles are considered, the depletion of the condensate has 
to be taken into account. The quantity 



n(0) = n c (0) + n(0) 



(9) 



is then a sum of the condensate density n c (0) and the 
density of noncondensed particles n(0). The latter term 



1357913579 

T (K) T (K) 

FIG. 6: The magnetization M(T) of TlCuCl 3 for 5.5 T 
< noH < 9 T applied along the &-axis (left) and along [201] 
(right), respectively. The critical temperature T c is marked 
by arrows. The solid lines correspond to fits to the data ac- 
cording to Eq. Q with S = 1 for (T < T c ). The dashed lines 
represent calculated M(T)-curves using Eq. and the same 
fitting parameters but with no=0. 

represents the number of triplons per Cu 2+ dimcr scat- 
tered out of the ground state due to the interactions be- 
tween the particles. It depends on the number of con- 
densed particles and can be expressed as @ 

1 { m*U o n c (0) \* 

n (°) = ^{— 2 — ) • (io) 

where we take the mass of a triplon m* w 2.6 x 10~ 29 
kg and the two-particle interaction potential C/o/^b ~ 
315 K from Ref.0. Replacing n(Q) in Eq. © with the 
expression in Eq. (|10[) we obtain 

1 fm*U o n c (0)\ ' , s 

n(0) =n c (0) + - ) ■ (ID 

From our fits according to Eq. ([7|) and with n(0) = 
^c«i + ri o we can now calculate the condensate density 
at zero temperature n c (0) for various magnetic fields us- 
ing Eq. flTjl. see Fig ETa). As one would expect from 
simple arguments [5j, [l9( ?t- c (0) increases with increasing 
magnetic field. It is essential to note that the number 
of triplons N C (Q) = n c (0) • Nd forming the condensate at 
T = is the same for both field directions only in the 
5=1 scenario for m up (T), see Fig. 13 and only in this 
scenario n c (0) extrapolates to zero at the correct critical 
field noH c w 5.5 T. These facts again strongly support 
our hypothesis that m up (T) is intrinsic with 5 = 1, and 
it confirms the interpretation of the magnetic field H act- 
ing as the chemical potential @ . 

In hqH = 6 T, right above the critical field H c , the per- 
centage of the condensed particles n c (0) with respect to 
the total density of triplons n(Q) is approximately 98% 
and slightly decreases with increasing magnetic field, see 
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(S — |), respectively. The error bars of Cs are of the order 
of the dot size. 

Fig. [5J This result is consistent with a similarly low non- 
condensed magnon density as calculated in Ref. [Til ], 
where n(0) increases from zero for H = H c to approxi- 
mately 7% of the total triplon density at T = in hqH = 
7 T. From the high percentage of condensed particles we 
can confirm that the triplons in TICUCI3 form a weakly 
interacting Bose gas [161 ] right above H c , and that the 
interaction increases with increasing particle density, i.e. 
with increasing the magnetic field H. 
Finally, we want to mention that the Curie-like contribu- 
tion Cs is small and essentially constant for H > H c in 
both the S = h and the 5=1 scenarios, see Fig. E^c). 
However, the corresponding data for 5=1 are more 
or less smooth continuations of the respective data for 
H < H c , in very contrast to the Cs data for 5 = | that 
show a discontinuity around H = H c , see Fig. [TJ^c). The 
comparably moderate variation of Cs with H over the 
whole considered range of magnetic fields for the 5=1 
scenario (covering both the normal phase and the BEC 
obeying an entirely different physics) may indicate that 
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FIG. 8: The condensate density n c (T = 0) with different 
scenarios for m up (left: S = 1, right: S — i). 
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FIG. 9: The triplon fraction n c (0)/n(0) forming the conden- 
sate in the S = 1 scenario for m UD . 



Cs is indeed a constant for each magnetic-field direction, 
and that our interpretation of a S = 1 paramagnetic 
background is correct. 



IV. CONCLUSIONS 

We have presented an analysis of magnetization 
M(T, H )-data of TICUCI3 and we calculated the density 
of condensed particles n c (0) at T = 0. The percentage of 
n c (0) with respect to the total density of triplons n(0) is 
approximately 98% in noH = 6 T and slightly decreases 
with increasing magnetic field. We demonstrated that 
this fraction is the same for both H || b and H || [201] 
if we assume the presence of a small number of intrin- 
sic S = 1 magnetic moments that are not part of the 
Bose-Einstein condensate of triplons even at the lowest 
temperatures. 




n H (T) 



FIG. 10: The critical triplon density n cr (a) and the corre- 
sponding typical temperatures T* « n cr Uo/kB < T c (b) in the 
S — 1 scenario for m up as functions of the applied magnetic 
field for H || b (circles and upward triangles) and H || [201] 
(squares and downward triangles), respectivley. 
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V. APPENDIX 

We want to emphasize that we do not interpret the 
exponent a in the power-law approach (Eq. Ij7]l) used for 
fitting the low-temperature magnetization data at high 
magnetic fields as a universal critical exponent. In this 
sense its physical meaning is not clear. 
In general, the normalfluid density in a dilute Bose gas in 
the condensed phase is proportional to T 4 at low enough 
temperatures (T <C T* ss nUo/kB), where n is the total 
particle density and Uq the interaction energy. 
In the case of TICUCI3 we can estimate T* by replacing 
n = n cr and Uo/ks ~ 315 K [i"9| . see Fig. [Till The values 
for a extracted from the fits, see Fig. [T^b) , vary around 
a ~ 4 with a decreasing tendency and increasing fitting 
error as H — > H c . This can be explained by the fact that 
for high magnetic fields T* is fairly close to T c , whereas 
for low magnetic fields the difference between the two 
characteristic temperatures increases, thereby restricting 
the validity of the T 4 - power law to very low tempera- 
tures that are not accessible in our experiment. Never- 
theless, it is clear that the increasing magnetization for 
T — > is related to the increasing fraction of the con- 
densed particles n c (T). The evaluation of n(0) from our 
phenomenological power law, see Eq. ([7]), and the calcu- 
lation of n c (0) using Eq. (fTTJ) gives, in any case, a reliable 
estimate of the intercept of n c (T) at T — 0, irrespective 
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of the correct functional form of n c (T). 
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